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Abstract. We present the results of a photometric sur- 
vey for variable stars in the nearby globular cluster M5. 
A 14.8 X 22.8 arcmin field centered on the cluster was 
monitored for a total of 37 hours with a CCD camera 
mounted on the 1.0-m Swope telescope. Five new vari- 
ables were identified: four blue stragglers which are SX 
Phe pulsating variables, and an eclipsing binary with an 
orbital period of 0.466 d. The eclipsing binary lies near 
the main-sequence turnoff on the cluster color-magnitude 
diagram. We do not confirm the variability of any of the 
10 eclipsing binaries identified in the cluster field by Reid 
(1996) and Yan & Reid (1996). The dwarf nova M5-V101 
exhibited two outbursts with full amplitudes of about 2.7 
mag during our observations. A V/B — V color-magnitude 
diagram of the surveyed field was obtained, and a possible 
extreme HB star located 2 mag below cutoff of the blue 
HB was identified. 

Key words: stars : binaries: eclipsing - blue stragglers - 
stars:variables:general - globular clusters: individual: M5 
- Hertzsprung-Russell and C-M diagrams 



1. Introduction 

NGC 5904=M5==C1516+022 is a rich globular clus- 
ter whose proximity and low reddening (r « 7.0 kpc, 
E{B - V) = 0.03; Harris 1996) make it an excellent 
target for detailed studies. We selected this cluster as one 
of the targets for an ongoing survey for eclipsing binaries 
in globular clusters (Kaluzny et al. 1997; Thompson 
et al. 1999). The main goal of this project is to use 
detached eclipsing binaries for the determination of ages 
and distances of globular clusters (Paczynski 1997). 

M5 is known to contain more than 120 RR Lyr vari- 
ables (Clement 1997 priv. comm.). CCD VI photometry 
for 49 of them was presented recently by Reid (1996) who 
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gives references to earlier papers dealing with RR Lyr stars 
in M5. As a side result of our survey we obtained F-band 
light curves for 65 RR Lyr stars from M5. These data will 
be presented in a separate contribution (Kaluzny et al. 
1999). 

Reid (1996) reported the discovery of 4 probable eclips- 
ing binaries among the cluster HB/RGB stars. In addition, 
Yan & Reid (1996) identified 6 short-period eclipsing bi- 
naries located at the cluster turn-off. Our data indicate 
that all of the eclipsing variables reported in Reid (1996) 
and Yan & Reid (1996) were spurious detections. 

Drissen & Shara (1998) used HST images to look for 
variable stars in the core of M5. They managed to identify 
one variable blue straggler, a possible eclipsing binary with 
a period P > 0.47 d. Several new RR Lyr stars were also 
identified in this study. 

M5 is one of only two galactic globular clusters known 
to harbor dwarf novae. Variable M5-V101 was first identi- 
fied as a dwarf nova by Oosterhoff (1941). A recent sum- 
mary of observational studies of that variable is given in 
Hakala et al. (1997). During our monitoring of M5 the 
variable exhibited two outbursts with full amplitudes ex- 
ceeding 2.7 mag in the V band. 



2. Observations and data reduction 

The time-series photometry of M5 on which this paper is 
based was obtained with the 1.0-m Swope telescope at Las 
Campanas Observatory. The data were collected during 
the interval 1997 May 9 - 1997 Aug 13. A condensed log 
of the observations is given in Table 1. 

A 2048 X 3150 pixel SITe CCD camera was used as 
the detector, with a field of view of 14.8 x 22.8 arcmin 
at a scale of 0.435 arcsec/pixel. Due to some technical 
problems the field of view of the camera was limited to 
14.8 x 14.8 arcmin during the period from 1997 May 9 
through 1997 May 19. During that run we monitored two 
overlapping subfields covering the northern and southern 
sections of the cluster. The total area of the monitored 
field was equal to 314 square arcmin. Starting in 1997 Jun 
7 we monitored just one 14.8 x 22.8 arcmin field centered 
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Table 1. Observing log 



Table 2. 

variables 



Equatorial coordinates for 5 newly identified 



Date 


Nv 


At 


UT 




hours 


May 9 


31 


5.5 


May 10 


34 


5.8 


May 13 


37 


5.6 


May 14 


37 


6.5 


May 15 


44 


5.8 


May 19 


24 


2.0 


Jun 7 


9 


2.0 


Jun 8 


12 


2.0 


Jun 17 


10 


1.5 


Aug 3 


1 


- 


Aug 4 


1 


- 


Aug 13 


2 


0.2 



on the cluster core. Monitoring was conducted in the V 
band with exposure times ranging from 300 to 500 sec- 
onds, depending on the seeing. The median value of the 
seeing for all of the images used for this survey was 1.52 
arcsec. Only two images in the B-band were taken dur- 
ing the 1997 season. On two nights in May/June 1998 we 
collected frames of the cluster field through the V and B 
filters using the same instrumentation as in 1997. These 
exposures were used to derive the color-magnitude dia- 
gram for the field monitored in 1997 (see Sect. 4). 

Instrumental photometry was extracted with 
DoPHOT (Schechter et al. 1993) used in the fixed- 
position mode. The stellar positions were derived from 
the reduction of "template" images, the best individual 
images obtained for a given sub-field. A total of 15839 
stars with V < 20.8 were checked for variability. The 
quality of the derived photometry is illustrated in Fig. 
1 in which we plot the rms deviation versus the average 
magnitude for stars from the northern sub-field. This plot 
includes 14847 stars with 13.5 < V < 20.8 whose light 
curves contain at least 75 points. Images of stars with 
V < 13.5 were generally over-exposed, and the resulting 
photometry was poor. To select potential variables we 
employed three methods as described in some detail in 
Kaluzny et al. (1996). Light curves showing possible 
periodic signals or smooth changes on time scales of days 
to weeks were selected for further examination. A total 
of 71 certain variables were identified: 65 RR Lyr stars, 
all of them already known, four previously unknown SX 
Phe stars, one previously unknown eclipsing binary, and 
a previously known cataclysmic variable (M5-V101). In 
Table 2 we list the equatorial coordinates for the 5 newly 
identified variables. Finding charts are shown in Fig. 2. 

2.1. Transformation to the standard system 

Observations of Landolt (1992) standards collected on sev- 
eral nights during the 1997 observing season were used 



ID 



RA(2000) 
h:ni:sec 



Dec(2000) 
deg:':" 



M5-NV1 15:18:22.76 02:02:49.3 

M5-NV2 15:18:11.40 02:00:40.6 

M5-NV3 15:18:40.04 01:58:41.6 

M5-NV4 15:18:41.02 01:56:29.6 

M5-NV5 15:18:39.10 02:11:41.0 



to establish transformations from the instrumental to the 
standard BV system. The following relations were ob- 
tained: 



y = ai + V- 0.010 x{B-V) + k.,xX 
b-v^a2 + 0.975 X {B-V) + kbyXX 

For the 1998 season we derived: 

v = ai+V - 0.006 x{B-V) + 0.126 x X 
b-v = a2+ 0.953 x {B - V) + 0.139 x X 



(1) 
(2) 

(3) 
(4) 



The extinction coefficients in Eqs. 3 and 4 correspond to 
the night of May 29/30, 1998, on which we obtained data 
used to construct the color-magnitude diagram for the 
field monitored for variables (see Sect. 4). While trans- 
forming the 1997 light curves of variables from the instru- 
mental V magnitudes to the standard V magnitudes we 
used a simplified relation: 



V = v + 2.435 



(5) 



This approach was required since no observations of the 
light curves in the B band were obtained. The constant 
offset was derived using photometry of non-variable stars 
with 14.5 < V < 17.0 and 0.1 < B - V < 0.75. The 
adopted procedure leads to systematic errors not exceed- 
ing 0.01 mag.[| 

2.2. Comparison with previous photometry 

The field covered by our survey overlaps the region stud- 
ied by Sandquist et al. (1996). Unfortunately, it does not 
include the field observed by Stetson and Harris (1988) 
and more recently by Johnson & Bolte (1998). Compari- 
son of our BV data with Sandquist et al. (1996) is shown 
in Fig. 3 (differences are given as "ours" - "theirs"). The 
mean differences for stars with V < 17.5 are 0.005 ±0.015 
and -0.021 ± 0.017, for V and B - V, respectively. A 
systematic trend oi AV residuals with increasing V mag- 
nitude can be noted in Fig. 3. A similar trend was noted 
by Johnson & Bolte (1998) who compared their VI data 
with Sandquist et al. (1996). 

^ The color term in Eq. 1 equals 0.010. All variables discussed 
in this paper have average colors in the range 0.20 < B — V < 
0.43. For RR Lyr variables changes of the B — V color seldom 
exceed 0.3 mag over the whole pulsational cycle. 



Table 3. Light curve parameters for newly identified vari- 
ables. The last column gives intensity averaged magnitude 
for SX Phe stars. 



ID 


P[day] 


B-V 


ymax 


V - 

y mm 


Av 


M5-NV1 


0.04213 


0.20 


17.42 


17.53 


17.48 


M5-NV2 


0.04130 


0.28 


16.75 


16.77 


16.76 


M5-NV3 


0.04166 


0.19 


16.99 


17.09 


17.04 


M5-NV4 


0.04745 


0.24 


16.89 


16.95 


16.92 


M5-NV5 


0.4663 


0.43 


18.91 


20.07 





3. Results for variables 

Our survey resulted in the detection of five new vari- 
ables: four SX Phe type stars and an eclipsing binary. 
Figure 4 shows the location of these objects on the cluster 
color-magnitude diagram (CMD). For the SX Phe stars 
the marked position corresponds to the intensity averaged 
magnitude and to the color measured near maximum light. 
For the eclipsing binary the marked position corresponds 
to the magnitude and color measured near the center of 
the secondary minimum at the orbital phase 0.52. The po- 
sition of the dwarf nova M5-V101 is also marked in Fig. 

4, and corresponds to the magnitude and color measured 
for that highly variable object on HJD=2450963.57 (May 
30, 1998 UT). 

The phased light curves of the four SX Phe variables 
are shown in Fig. 5 and the light curve of the eclipsing 
binary is shown in Fig. 6. Some basic photometric param- 
eters of the newly identified variables are listed in Table 3. 
The periods of variability were derived using an algorithm 
based on the analysis oj variance sta,t\si\c (Schwarzenberg- 
Czerny 1989). 

3.1. SX Phe stars 

The SX Phe variables are located in the blue straggler 
region of the cluster CMD. The observed properties of 
these variables allow their classification as SX Phe stars. 
Fig. 7 shows the positions of the variables in a period vs. 
absolute magnitude diagram. The standard relations for 
SX Phe stars pulsating in the fundamental mode and in 
the first overtone are also shown (McNamara 1997, sec also 
McNamara 1995). The plotted relation for stars pulsating 
in the fundamental mode is 



Mv 



-3.725 X LogPo - 1-933 



(6) 



and the same relation is used for stars pulsating in the 
first overtone assuming Pi/Pq — 0.778. We also plot the 
P — L — [Fe/H] relations for fundamental and first over- 
tone pulsators from Nemec et al. (1994), assuming metal- 
licity of [Fe/H] = —1.4. The absolute magnitudes for M5 
variables were calculated assuming a distance modulus of 
(m - M)v = 14.33 (Harris 1996). It is clear that the only 
variable in Fig. 7 close to the standard relation for the fun- 
damental mode pulsators is M5-NV1. It is likely that the 



3 remaining SX Phe stars are first overtone pulsators. The 
small amplitudes and the symmetry of their light curves 
support such a hypothesis, though as McNamara (1997) 
points out, classification of a star as a first overtone pul- 
sator based only on these criteria is questionable. How- 
ever, even assuming that the 3 brightest of the 4 iden- 
tified SX Phe stars are first overtone pulsators, we still 
face the problem of an apparent over-luminosity of all 4 
variables most especially for the Nemec et al. (1994) rela- 
tions. One possible solution to this problem is adoption of 
a lower distance modulus for the cluster, {■m—M)v ~ 14.1, 
a value entirely discrepant with that adopted by Harris 
(1996) for this well studied cluster. This would also imply 
a very low luminosity for the cluster RRab variables of 
My ~ +1.00 (the average V magnitude for RRab stars 
in M5 is 15.062 (Kaluzny et al. 1999)). Another possibil- 
ity is that discovered SX Phe stars are not members of 
M5, it although M5 is at a galactic latitude of 5 = 46.8 
deg, and this expalnatgion implies un unreasonably high 
field surface density for SX Phe stars. One way to resolve 
this puzzle would be to measure the radial velocities for 
these variables to check their membership status. A more 
likely explanation is that the absolute magnitude - period 
relations are poorly known for stars of such short peri- 
ods. Indeed, McNamara's calibration based on Hipparcos 
measurements has only five stars with periods below log 
P = —1.1, and the scatter in this region of the calibration 
is significant. It is hoped that our survey for variable stars 
in other globular clusters will lead to the discovery of sig- 
nificant populations of SX Phe stars, leading to a more 
secure calibration. 



3.2. Eclipsing binary stars 

Variable M5-NV5 is potentially an important object as 
it may prove to be a detached eclipsing binary composed 
of two upper-main sequence stars from M5. As such this 
binary could be used for measuring the distance to the 
cluster. In addition, the determination of the masses, radii 
and luminosities of the components of this system is of 
great interest. However, the present data are not sufficient 
for a precise determination of parameters of M5-NV5. 

To get some hints about the membership status and 
configuration of M5-NV5, we analyzed the light curve pre- 
sented in Fig. 6. The light curve was solved using the 
MINGA package (Plewa 1988) which uses the Wilson- 
Devinney (W-D) code (Wilson & Devinney 1971, Wilson 
& Sofia 1976 ) to generate synthetic light curves. The code 
was used in a mode in which the basic parameters deter- 
mining the shape of the synthesized light curve are the 
mass ratio q = M2/M1, the potentials fli and ^2 (for 
a given q these determine the relative radii ri and r2), 
the relative luminosities Liy and L2V, and the inclina- 
tion of the orbit i. During the analysis we fixed the tem- 
perature of the primary component at Ti = 6300 K. Ac- 
cording to the derived solution the secondary eclipse is 
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Table 4. Light curve solutions for M5-NV5. 



Q 


i 


L\v / L2V 


ri 


r2 


0.50 


89.2 


3.68 


0.32 


0.27 


0.60 


88.9 


3.78 


0.32 


0.27 


0.70 


87.2 


3.82 


0.31 


0.27 


0.80 


87.8 


4.07 


0.32 


0.27 


0.90 


85.9 


3.45 


0.30 


0.28 



total or very close to total. Hence at an orbital phase of 
0.52, at which we measured B — V, only the primary com- 
ponent is visible. Adopting E{B — V) = 0.03 we obtain 
{B — V)i ~ 0.40 which in turn implies Ti = 6300 for an 
adopted [Fe/H] = -1.1 (Alonso et al. 1996; Harris 1996). 
As we have no information about the mass ratio of the 
system we calculated a grid of light curve solutions for a 
set of values of q. The results are summarized in Table 4. 

All of the solutions in Table 4 imply that the system 
is detached with the secondary eclipse being total or very 
close to total. Based on a x^ criterion, the best solution 
was obtained for q — 0.70, but solutions derived for other 
values of q are only marginally worse. This is not surpris- 
ing considering the limited accuracy of the photometry 
and the detached configuration of the binary. To illus- 
trate the quality of the solutions we plot in Fig. 8 the 
observed light curve with the synthetic light curve corre- 
sponding to q = 0.70. Despite our inability to derive the 
mass ratio of the system it is safe to state that M5-NV5 
is a detached binary with an inclination -i > 86 deg and 
that the luminosity ratio Lw / L2V > 3.4. The position of 
the primary component of M5-NV5 on the cluster CMD 
(see Fig. 4) is consistent with the hypothesis that it is a 
main-sequence star belonging to the cluster. Finally, we 
note that the derived values of Lw / L2V and ri/r2 favor 
solutions with 0.7 < q < 0.8 (if we assume that both com- 
ponents are main sequence stars). If the system contains 
two cluster main sequence stars with masses toi « 0.7Mq 
and TO2 ~ 0.5Mq then the expected semi-amplitudes of 
their radial velocity curves are Ki ~ 121 and K2 = 170 
km/sec. In this case the velocity curve of the binary would 
be relatively easily measured with a medium resolution 
spectrograph on a large telescope. 

In a paper devoted mainly to a photometric analysis 
of RR Lyr stars in M5, Reid (1996) reported the identifi- 
cation of 4 eclipsing binaries - 2 located on the horizontal 
branch and 2 located on the giant branch. Subsequently 
Yan & Reid (1996) used the same observational material 
to identify 6 short-period eclipsing binaries located on or 
near the main-sequence of the cluster. None of these 10 
candidate variable stars was selected from our data base 
during our search for variable objects. Our light curves 
for these stars do not show any sign of variability of the 
type reported by Reid (1996) and Yan & Reid (1996). At 
our request Dr. Reid checked his CCD frames of M5. It 
turned out that the variability of the 4 candidate eclips- 



ing binaries reported by Reid (1996) was due to occasional 
passages of the stellar images through defective columns 
of the detector. Our light curves for the 6 candidate vari- 
ables reported by Yan & Reid (1996) show scatter ranging 
from 0.04 mag for V6 and V5 to 0.08 mag for V2. None of 
these light curves show any modulation with the periods 
or at the amplitudes reported by Yan & Reid (1996), nor 
with any other periods shorter than the time-base of our 
observations. This is demonstrated in Fig. 9 in which we 
show our light curves for stars VI, V2 and V3 from Yan & 
Reid (1996) phased with the periods listed in that paper. 

3.3. The cataclysmic variable M5-V101 

M5-V101 is one of two cataclysmic variables which were 
known in globular clusters in the pre-HST era. It was dis- 
covered by Oosterhoff (1941), who suggested that it is a 
dwarf nova. That suggestion was later confirmed by spec- 
troscopic as well as photometric observations (Margon et 
al. 1981; Shara et al. 1987). Detection of the X-ray coun- 
terpart of M5-V101 was reported by Hakala et al. (1997). 
In Fig. 10 we present the V^-band light curve of M5- 
VlOl based on our observations. The variable displayed 
two outbursts with full amplitudes of about 2.7 mag dur- 
ing our observations in 1997. On the night of May 30 
1998 we observed M5-V101 at V = 20.27 ± 0.04 and 
B -V = 0.10 ± 0.08 (HJD = 245 0963.57). The occur- 
rence of at least two outbursts on a time scale of 100 days 
in 1997 as weU as the detection by Shara et al. (1987) of a 
single outburst during a two- week long monitoring period 
indicates a relatively short duty cycle for M5-V101. Exam- 
ination of the nightly light curves from the 1997 season re- 
veals occasional variability with a full amplitude of about 
0.4 mag on a characteristic time scale of about 3.4 hours. 
Spectroscopic data presented by Naylor et al. (1989) sug- 
gest that the orbital period of M5-V101 is longer than 3.5 
hours. 

4. The color-magnitude diagram 

An extensive and detailed photometric study of M5 based 
on CCD BVI data was published by Sandquist et al. 
(1996). Even more recently Johnson & Bolte (1998) re- 
ported VI data for a field in the outer part of the cluster. 
As a by product of our program we obtained a medium- 
deep VI B - V CMD for the 14.8 x 22.8 arcmin field cen- 
tered on the cluster, shown in Fig. 4. This CMD is based 
on the following set of images: i? - 60 s, 150 s, 260 s; 
V - 60 s, 4 X 150 s. The observations were obtained on 
the nights of 1998 May 30 and June 3. The seeing for 
these frames ranged from 1.13 to 1.35 arcsec. Photome- 
try was extracted using Daophot/AUstar package (Stetson 
1987). Measurements of stars with relatively poor quality 
(large sigma for a given magnitude and/or a large value of 
Stetson's CHI parameter) were rejected from the analysis. 
Known variables listed in Clement (1997, priv. comm.) or 



in Sandquist ct al. (1996) are not plotted in Fig. 4. The 
CMD shows several features discussed in some detail by 
Sandquist et al. (1996). In particular we note the pres- 
ence of a star lying to the blue of the HB at V — 15.85 
and B — V = —0.27. This star is located at angular dis- 
tance r = 10.1 arcmin from the projected cluster center. 
This does not exclude its cluster membership but clearly 
additional data are needed to clarify this point. The equa- 
torial coordinates of this star are 02000 = 15 : 18 : 32.7 
and (52000 = 01 : 54 : 52.8. The CMD also includes a 
faint blue star with V = 18.836 and B - V ^ -0.306. 
It is located about 2 magnitudes below the cutoff of the 
blue horizontal branch of the cluster. This star is marked 
with the triangle in Fig. 4. Its equatorial coordinates are: 
a2ooo = 15'' : 18™ : 14.5^ and (52ooo = 02° : 04' : 09.8". 
It is a good candidate for an extreme horizontal branch 
star. The light curve of this star shows no evidence for 
variability with a full amplitude exceeding 0.02 mag. 

5. Conclusions 

The most interesting result of this paper is the identifica- 
tion of a detached eclipsing binary which is most probably 
composed of two upper main sequence stars from M5. Pho- 
tometric and spectroscopic follow up observations of this 
star are planned by our group. Our data indicate that 10 
candidate eclipsing binaries reported by Reid (1996) and 
Yan & Reid (1996) do not show any periodic variability 
on time scales shorter than about a week. Four SX Phe 
stars were identified, all located among the blue stragglers 
on the cluster CMD. Assuming cluster membership, their 
absolute magnitudes My are 0.2-0.3 mag too bright in 
comparison with the standard relations obtained by Mc- 
Namara (1997, 1995). Two outbursts of the dwarf nova 
M5-V101 were detected during the period May- Aug 1997. 
A faint blue object, which we propose is an extreme hori- 
zontal branch star from M5, was identified. 
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Fig. 1. Standard deviation versus average V magnitude 
for 14847 stars from the north field. Only stars measured 
on at least 75 out of the 138 best frames are marked. 
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Fig. 2. Finder charts for newly identified variables. Each 
chart is 44 arcsec on a side with north up and east to the 
left. 
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Fig. 3. Comparison of our photometry with Sandquist et 
ah (1996). 
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Fig. 4. A V/B - V CMD for M5 based on our data. Filled 
triangles correspond to SX Phe stars. The position of the 
eclipsing binary is marked with a square. Known RR Lyr 
stars are not plotted. The position of M5-V101 is marked 
with an asterisk 
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Fig. 5. Phased light curves for 4 SX Phe stars from M5. 
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Fig. 6. Phased light curve for the echpsing binary star 

M5-NV5. 
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Fig. 7. The period - luminosity diagram for SX Phe stars 
from the field of M5. Standard relations for fundamental 
(solid line) and first overtone pulsators (dashed line) are 
marked: M - McNamara (1997); N - Nemec et al. (1994). 
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Fig. 8. The observed light curve of M5-NV5 overplotted 
with a synthetic hght curve corresponding to an assumed 
mass ratio oi q = 0.70. 
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Fig. 9. Observed light curves for stars VI (top), V2 and 
V3 (bottom) from Yan & Reid (1996) phased with periods 
listed in that paper. 
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Fig. 10. The V band light curve of the dwarf nova M5- 
VlOl. 
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Fig. 11. A V/B-V CMD for the monitored field. The 
trangle corresponds to a possible cluster extreme horizon- 
tal branch star. 



